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♦♦♦ Leading-edge load 
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P What the technology does 


Problem Statement 

□ Improving fuel efficiency for an aircraft 

❖ Reducing weight or drag 

> Similar effect on fuel savings 

❖ Multidisciplinary design optimization [design phase) or active control [during 
flight) 

□ Real-time measurement of deflection, slope, and loads in flight are a valuable tool. 



Complete degrees of freedom 

Deflection 


Slope [angle) 


Wires for Strain Gage Wire for FOSS 


□ 


□ 


Wing deflection and slope [complete degrees of freedom) are essential quantities 
for load computations during flight. 

Loads can be computed from the following governing equations of motion. 

M q{t) + G q(t) + K q(t) = F a {Mach,q{t )) 

> Internal Loads: using finite element structure model 

S Mq(t), G q(t), K q(t) : Inertia, damping, and elastic loads 

> External Load: using unsteady aerodynamic model 
S F a : Aerodynamic load 

Traditionally, strain over the wing are measured using strain gages. 

❖ Cabling would create weight and space limitation problems. 

❖ A new innovation is needed. Fiber optic strain sensor [FOSS) is an ideal 
choice for aerospace applications. 



FOSS „ . 

Strain Gage 


Structural Dynamics Group 


Wing deflection & slope will be computed from measured strain. 
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Previous technologies 


Beam theory; Sectional bending moment and shear loads 


□ Liu, T., Barrows, D. A., Burner, A. W., and Rhew, R. D., "Determining Aerodynamic Loads Based on Optical Deformation Measurements,” AIAA 2001-0560 

❖ NASA LRC; Application is limited for "beam". 

□ Shkarayev, S., Krashantisa, R., and Tessler, A., "An Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and 
Structural Response of Aerospace Vehicles,” Proceedings of Third International Workshop on Structural Health Monitoring, 2001 

❖ University of Arizona and NASA LRC; using an inverse interpolation formulation. 

□ Kang, L.-H., Kim, D.-K., and Han, J.-H., "Estimation of Dynamic Structural Displacements using fiber Bragg grating strain sensors,” 2007 

❖ KAIST; displacement-strain-transformation (DST) matrix. Use strain mode shape. Application was based on beam structure. 

□ Igawa, H. et al., "Measurement of Distributed Strain and Load Identification Using 1500 mm Gauge Length FBG and Optical Frequency Domain 
Reflectometry,” 20th International Conference on Optical Fibre Sensors, 2009 

❖ JAXA; using inverse analysis. "Beam" application only. 

□ Ko, W. and Richards, L., "Method for real-time structure shape-sensing,” US Patent #7520176B1, April 21, 2009 

❖ NASA AFRC; closed-form equations (based on beam theory) 

□ Richards, L. and Ko, W. , "Process for using surface strain measurements to obtain operational loads for complex structures,” US Patent #7715994, May 

11,2010 

❖ NASA AFRC; "sectional" bending moment and shear force along the "beam". 

□ Moore, J.P., "Method and Apparatus for Shape and End Position Determination using an Optical Fiber,” U.S. Patent No. 7813599, issued October 12, 2010 

❖ NASA LRC; curve-fitting 

□ Park, Y.-L. et al., "Real-Time Estimation of Three-Dimensional Needle Shape and Deflection for MRI-Guided Interventions,” IEEE/ASME Transactions on 
Mechatronics, Vol. 15, No. 6, 2010, pp. 906-915 

❖ Harvard University, Stanford University, and Howard Hughes Medical Institute; Uses beam theory. 

□ Carpenter, T.J. and Albertani, R., "Aerodynamic Load Estimation: Pressure Distribution from Virtual Strain Sensors for a Pliant Membrane Wing,” AIAA 
2013-1917 

❖ Oregon State University; Aerodynamic loads are estimated from measured strain using virtual strain sensor technique. 
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Previous technologies are applied to a beam structure. 
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Technical features of new technology 


Proposed solutions: 

□ The new method for obtaining the deflection over a flexible full 
3D aircraft structure is based on the following two steps. 

❖ First Step: Compute wing deflection along fibers using measure 
strain data 

> Wing deflection will be computed along the fiber optic sensor line. 

> Strains at selected locations will be "fitted”. 

> These fitted strain will be integrated twice to have deflection 
information. (Relative deflection w.r.t. the reference point) 

> This is a finite element model independent method. Drag and 

❖ Second Step: Compute wing slope and deflection of entire structures lift 

> Slope computation will be based on a finite element model 
dependent technique. 

> Wing deflection and slope will be computed at all the finite element 
grid points. 
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A new two-step theory is investigated for predicting the deflection and slope of an entire structure using strain measurements at discrete locations. 
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Piecewise least squares curve fit boundaries 


Technical features of new technology (continued) 

□ First Step 

❖ Use piecewise least-squares method to minimize noise in the 
measured strain data (strain/offset) 

❖ Obtain cubic spline (Akima spline) function using re-generated strain 
data points: 

d 2 S 

d . ? = “ e(s)/c(s) 

❖ Integrate fitted spline function to get slope data: 

❖ Obtain cubic spline (Akima spline) function using computed slope data 

❖ Integrate fitted spline function to get deflection data: 

S(s) 


-.007 



o : raw data 

— : direct curve fit 
: curve fit after piecewise LS 

— 
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Along the fiber direction, in. 
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A measured strain is fitted using a piecewise least-squares curve fitting method together with the cubic spline technique. 
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Technical features of new technology (continued) 

Second Step: Based on General Transformation 

❖ For all model reduction/expansion techniques, there is a relationship between the master (measured or tested) degrees of freedom and the 
slave (deleted or omitted) degrees of freedom which can be written in general terms as 


(«} = {’"} = mra*}: 


[q}= general displacement vector 

to); 


S Where, an eigen-matrix is defined as j = 


0 


si 


{rj}= orthogonal displacement vector 


Transformation matrix [T] can be one of the followings: 

S Guyan (or static) condensation, dynamic condensation, improved reduced system (IRS), or system equivalent reduction expansion 
process (SEREP) 


□ Expansion of displacement using SEREP: kinds of least-squares method; most accurate reduction-expansion technique 
❖ {qr M }: master DOF; deflection along the fiber "computed from the first step" 



Computed deflection along the fibers are combined with a finite element model of the structure in order to interpolate and extrapolate 
the deflection and slope of the entire structure through the use of the System Equivalent Reduction and Expansion Process. 



Computational Validation 



Cantilevered rectangular wing model 



Cantilevered Rectangular Wing Model 


□ Wind tunnel test wing (thickness = 0.065 in.) 

❖ Uniform lg load 


Uniform lg loading 


❖ Wing tip torsion (1 lbf at leading-edge and -1 lbf at trailing-edge of wing tip section) Is 


♦♦♦ Aerodynamic load under 1° angle of attack at Mach 0.715 

□ MSC/NASTRAN 

❖ Compute strain 

❖ Compute deflection (target) 

□ ZAERO 

❖ Compute aerodynamic load 

□ Two-step approach 

❖ Compute deflection from computed strain 

❖ Compare computed deflection with respect to target value 


Strain (y component) ^ 



Wing tip torsion 

(1 lbf at leading-edge and -1 lbf at trailing-edge of wing tip section) 
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Aerodynamic load under 1° angle of attack at Mach 0.715 


Plate 
21 elements 
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Cantilevered Rectangular Wing Model: Uniform lg 


□ Uniform lg load 



Wing deflection over FE model 
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Along the fiber direction, in. 


b) Wing deflection (step 2 results with 10 modes) 



a) Curvature distribution (step 1 results) 
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d) Wing span-wise slope (step 2 results with 10 modes) 

Chan-gi Pak-10 



Cantilevered Rectangular Wing Model: Wing tip torsion 


□ Wing tip torsion 
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Positive 
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Markers: NASTRAN Result Lines: Two-step Approach 
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(d) Step 2 results with 50 Modes 
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Along the Fiber Direction (inch) 


Markers: NASTRAN Result Lines: Two-step Approach 
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(a) Strain distribution (step 1 results) 
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(d) Step 2 results with 50 Modes 


(d) Step 2 results with 50 Modes 
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Curvature, /in. 



Cantilevered Rectangular Wing Model: Aerodynamic load 


□ Aerodynamic load under 1° angle of attack at Mach 0.715 



Markers: NASTRAN Result 





a) Curvature distribution (step 1 results) 
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(b) Splined load vectors 
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Slope (Roll Direction) X Deflection (inch) 


y Cantilevered Rectangular Wing Model: Aerodynamic load (continued) 
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Markers: NASTRAN Result Lines: Two-step Approach 


Markers: NASTRAN Result Lines: Two-step Approach 
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X deflection (step 2 results with 10 modes) 
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Markers: NASTRAN Result 
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Y deflection (step 2 results with 10 modes) 


Markers: NASTRAN Result Lines: Two-step Approach 
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Z deflection (step 2 results with 10 modes) 
Markers: NASTRAN Result Lines: Two-step Approach 
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Along the Fiber Direction (inch) 

Wing slope in roll direction (step 2 results with 10 modes) 
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Along the Fiber Direction (inch) 

Wing slope in pitch direction (step 2 results with 10 modes) 
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Along the Fiber Direction (inch) 

Wing slope in yaw direction (step 2 results with 10 modes) 
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p Cantilevered Rectangular Wing Model: Aerodynamic load (continued) 


□ Aerodynamic load under 1° angle of attack at Mach 0.715 [continued) 

□ Wing tip deflections {q s } : slave DOF 

Deflection in X. Y. & Z direction 


Boundary 


Fiber number 

Target (inch) 

Computed X deflectidnlinchl 

Relative error (%) 

Step 1 

Step /with 1 (’hpodes 

Step 1 

Step 2 with 10 modes ^ 

i 

0.0006927 

N/A 

/ 0.0006703\ 

N/A 

d -3.2 

5 

0.0006975 

N/A 

/ 0.0006776 \ 

N/A 

-2.8 

9 

0.0007029 

N/A 

/ 0.0006936 \ 

N/A 

-1.3 

13 

0.0007047 

N/A 


0.0007081 \ 

N/A 

0.48 

17 

0.0007035 

N/A 


0.0007144 


N/A 

1.6 

21 

0.0007023 

N/A 


0.0007143 


N/A 

1.7 

Fiber number 

Target (inch) 

Computed 

{ deflection (inch) 


Relative error (%) 

Step 1 


Step 2 with 10 mode; 


Step 1 

Step 2 with 10 modes 

i 

-0.003280 

N/A 


-0.003282 


N/A 

0.08 

5 

-0.003272 

N/A 


-0.003278 j 


N/A 

0.18 

9 

-0.003260 

N/A 

\ -0.003268 / 

N/A 

0.26 

13 

-0.003248 

N/A 

\ -0.003257 / 

N/A 

d 0.28 3 

17 

-0.003239 

N/A 

\ -0.003246 / 

N/A 

0.23 

21 

-0.003235 

N/A 

\0.00324/ 

N/A 

0.20 

Fiber number 

Target (inch) 

Computed Z deflectjs^mch) 

Relative error (%) 

/Step T\ 

StepyZ^with luSnodes 

Step 1 

Step 2 with 10 modes 

i 

0.9280 

/ 0.9269 \ 

j 0.9275 \ < 

C3).12 ^ 

> -0.06 

5 

0.9085 / 0.9091 \ 

I 0.9087 \ 

0.07 

0.02 

9 

0.8889 0.8894 

1 0.8895 1 

0.06 

0.08 

13 

0.8691 l 0.8696 / 

\ 0.8699 / 

0.05 

d 0.09 3 

17 

0.8493 

\ 0.8497 / 

V 0.8498 y 

0.05 

0.06 

21 

0.8296 

A^830/ 

NCL822X 

0.04 

0.01 



Undeformed 

shape 


Wing deflection over FE model 
fstep 2 results with 10 modesl 


Smoothing effect 
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Input to Step 2 master DOF {q M }: smoothed master DOF 
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p Cantilevered Rectangular Wing Model: Aerodynamic load (continued) 


□ Aerodynamic load under 1° angle of attack at Mach 0.715 [continued) 

□ Wing tip slopes {qs} . slave D0F 

Slope in roll, pitch. & yaw direction 


Not usable for Step 2 computations 


Y deflection 


Fiber number 

Target 

Computed \)11 slope 

/ Relative error (%) 

Step 1* 

Step 2ydTiri>Q modes 


Step 2 with 1 0 mod 

i 

0.10090 

0.1010 

yC.lOlOoX 

f 1.12^ 

/t).08' v \i^ 

5 

0.10070 

0.1016 

/ 0.10090 \ 

/ 0.96 \ 

/ 0.18 \ 

9 

0.10030 

0.1012 

/ 0.10060 \ 

0.94 

f 0.26 \ 

13 

0.09993 

0.1009 

/ 0.10020 \ 

0.93 1 

ill 0.28 Zp 

17 

0.09966 

0.1006 


0.09989 \ 


V 0.92 / 

\ 0.23 j 

21 

0.09954 

0.1004 


0.09974 


\0.88 / 

\^.20/ 

Fiber number 

Target 

Comput 

:d pitch slope 


Relative error (%) 

Step 1 


itep 2 with 10 mod( 

s 

Step 1 

Step 2 with 10 modes 

i 

0.02131 

N/A 


0.02063 


N/A 

CZ-T2>0 

5 

0.02146 

N/A 


0.02085 


N/A 

7 -2.9 T 

9 

0.02163 

N/A 


0.02134 


N/A 

1 ^ 1 

13 

0.02168 

N/A 


0.02179 


N/A 


17 

0.02165 

N/A 


0.02198 


N/A 

Y 1.5 / 

21 

0.02161 

N/A 


0.02198 


N/A 

\T77^n 

Fiber number 

Target 

Compu 

ed yaw slope 


Absolute error . 

Step 1 

< 

tep 2 with 10 mod 

:s 

Step 1 

Step 2 with 10 modes 

i 

2.2e-31 

N/A 


7.9e-18 / 


N/A 

0.0000 

5 

1.9e-31 

N/A 

\ 6.6e-18 / 

N/A 

0.0000 

9 

1.7e-31 

N/A 

\ 6. 1 e-1 8 / 

N/A 

0.0000 

13 

1.7e-31 

N/A 

\ 5.7e-18 / 

N/A 

0.0000 

17 

1.4e-31 

N/A 

\ 5.9e-18/ 

N/A 

0.0000 

21 

1.2e-31 

N/A 

V.3e-iX 

N/A 

0.0000 


Relative error (%) 

Step 1 

Sffep 2 with 10 modes 

N/A 

y'txoK 

N/A 

/ 0.18 \ 

N/A 

( 0.26 \ 

N/A 

\ 0.28 1 

N/A 

\ 0.23 / 

N/A 

V0.2D/ 

Relative error (%) 

Step 1 

Step 2 with 10 modes 

N/A 


N/A 

/ -2.8 \ 

N/A 

[ U3 | 

N/A 

“T 0.48 1 

N/A 

Y 16 / 

N/A 



X deflection 


: Roll slope without effect of X and Y deflections (These slopes are not used during step 2 computation.) 
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Experimental T esting 



Swept test plate 


V 


Swept Test Plate 

□ Tested at NASA AFRC 



a) Leading-edge load b) Uniform load 
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Curvature, / in. 


Swept Test Plate (continued) 

□ Averaging the curvatures calculated by using each fiber individually eliminates the effect of the axial load. 

□ This computation is performed after curve-fitting each set of data individually to minimize noise. 




(a) Leading-edge load 


(b) Uniform load 
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Average deflection, in. Average deflection, in. 



Swept Test Plate (continued) 


□ Under leading-edge load 
o.o 



O : Photogrammetry data 
— : Computed deflection 
- : Bakalyar/Jutte deflection 

— 
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Along the fiber direction, in. 

a) Leading-edge fiber (step 1 result) 
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O : Photogrammetry data 
— : Computed deflection 
-- - : Bakalyar/Jutte deflection 
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Along the fiber direction, in. 
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b) Leading-edge fiber (step 2 result with 10 modes) 


d) Mid-chord fiber (step 2 result with 10 modes) 


f) Trailing-edge fiber (step 2 result with 10 modes) 
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Average deflection, in. Average deflection, in. 


Swept Test Plate (continued) 


□ Under uniform load 


o.o 
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Step 1 results 


O : Photogrammetry data 
— : Computed deflection 
: Bakalyar/Jutte deflection 


- 


- 


Leading-edge 


10 20 30 

Along the fiber direction, in. 

a) Leading-edge fiber (step 1 result) 
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Mid-chord 


O : Photogrammetry data 
— : Computed deflection 
: Bakalyar/Jutte deflection 


10 20 30 

Along the fiber direction, in. 

c) Mid-chord fiber (step 1 result) 
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larger error 


O : Photogrammetry data 
— : Computed deflection 
: Bakalyar/Jutte deflection 






Trailing-edge 
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Along the fiber direction, in. 

e) Trailing-edge fiber (step 1 result) 
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Along the fiber direction, in. 


f) Trailing-edge fiber (step 2 result with 10 modes) 
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Swept Test Plate (continued) 


□ Deformed wing shape (step 2 results with 10 modes) 



Measured 

(inch) 

Computed (inc 

h) 

Relative error (%) 

Bakalyar 
and Jutte* 

Step 1 

Step 2 
with 10 
modes 

Bakalyar 
and Jutte 

Step 1 

Step 2 
with 10 
modes 


^eading-edge load / — \ /~\ 

Leading edge fiber 

-4.525 

-4.500 

-4.542 

-4.569 

/ -0.55 \ 

f 0.38\ 

0.97 

Middle fiber 

-4.912 

-4.952 

-4.880 

-4.843 

/ 0.81 \ 

/ -0.65 \ 

-1.40 

Trailing edge fiber 

-5.300 

-5.067 

-5.091 

-5.097 


-4.40 

-3.90 

-3.sn^> 

Uniform load 




Leading edge fiber 

-6.541 

-6.546 

-6.630 

-6.684 


0.08 / 

1.40 

2.20 

Middle fiber 

- 7.256 

-7.408 

-7.313 

- 7.238 

\ 210 / 

/C 0.79 7 

-CL25-N 

Trailing edge fiber 

-7.971 

-7.667 

- 7.750 

- 7.763 

\-3.80 / 

V\_-2.80/ 

-2.60^ 


*. 


extrapolated result Boundary 


Boundary 



Undeformed 
shape 

Deformed 
shape 

r a) Leading-edge load. 
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Undeformed 
shape 



In general larger error than current approach 


Deformed 

shape 


b) Uniform load. 


Smoothing effect 



Stress concentration 


Boundary 
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Conclusions 

□ The two-step approach for computing all the degrees of freedom in a structural FE model from 

measured strain along the FOSS is successfully applied to a cantilevered rectangular wing model and a 
test plate. 

❖ The first experiment investigates the accuracy of the theory by applying it to a cantilevered 
rectangular wing model analyzed using the MSC/NASTRAN and ZAERO codes. 

> lg uniform load case 

> Wing tip torsion load 

> Aerodynamic loading 

S All six computed DOFs have excellent matching with target values. 

❖ The second experiment applies the theory to experimental data collected from a test plate 
fabricated and tested at the NASA AFRC. 

> The deflections calculated from the experimental model are extremely accurate. 
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Average Deflection (inch) 


Boundary 
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